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Abstract. It is shown that the slow ghtchcs in the spin 
rate of the pulsar B1822-09 can be explained by the recon- 
struction of the neutron star shape, which is not matched 
with the star rotation axis. Owing to the evolution of the 
inclination angle, i.e. the angle between the rotation axis 
and the axis of the magnetic dipole, under the action of 
the braking torque, there appears the disagreement be- 
tween the rotation axis and the symmetry axis. After the 
angle between the axis of symmetry and the axis of the 
rotation achieves the maximum value of a ~ 2 ■ 10~^ the 
shape of the neutron star becomes matching with the ro- 
tation axis. Such reconstruction is observed as the slow 
glitch. 



1. Introduction 

All radio pulsars are known to retard theirs rotation. Pe- 
riods of the rotation P are gradually increasing with the 
rate of the order of P « 10~^^s/s. As a result during the 
time P/P « lO'^ years the rotation period becomes twice 
larger. At the same time the inclination angle x ~ the an- 
gle between the rotation axis and the axis of the magnetic 
dipole - must changes also. The characteristic time of the 
evolution of the x is the same as for the rotation period 
P. The change of the inclination angle is of « 1 during 
the evolution time. Change of the value of x is due to the 
torque acting on the star is not parallel to the rotation 
axis. This torque not only retards the rotation but also 
turns the rotation axis in space. Increasing or decreasing of 
the X depends on the mechanism of a neutron star energy 
losses. The angle x approaches to the value when the en- 
ergy losses are minimal. For the magnetodipole losses, if a 
pulsar radiates the electromagnetic magnetodipole waves, 
the angle x goes to the zero value when the radiation is 
absent. Opposite, for the current losses, if a pulsar spends 
its energy for the creation and acceleration of plasma in a 
magnetosphere, and a star's spin down is due to the elec- 
tric currents flowing in a magnetosphere and closing on a 
star's surface, the angle x approaches to the value of 90°. 
In this case the electric current is minimal. 



However, it is clear that to notice during the obser- 
vation time w 10-20 years the change of the inclination 
on the value of 10~^rad « 0."2 is impossible. Thus, there 
exist only indirect conclusions about inclination angle evo- 
lution. They are based on the statistic of the distribution 
of pulsars over the value of x- But because we do not know 
this distribution at the birth time, and due to the selec- 
tion effects we can't now do the definite conclusion about 
inclination angle evolution. 

Our work demonstrates the possibility to watch the 
changes of x on the value of 2' during the time « lOyears 
though slow ghtches observed from the pulsar B1822-09. 

2. Observed properties of slow glitches from the 
pulsar B1822-09 

The pulsar B1822-09 has the period P = 0.769s = 
1/P = 1.3s-i), the period derivative P = 52.36 x 10"^^. It 
is the young pulsar with the age of P/2P « 2.3 x lO^years. 
The pulsar observations during the time of 21 years at 
Puschino Radio Observatory show two class of the period 
variations. They are pulsar noise and glitches. The pulsar 
B1822-09 reveals the new kind of glitches unknown before 
- slow glitches (Shabanova 2005). 

Figure 1 shows the frequency first derivative v and 
frequency residuals Aiy for PSR B1822-09 from 1985 to 
2006. The local fits were performed over intervals of 200 
days. The figure 1(b) presents the frequency residuals with 
respect to the simple spin down model based on the initial 
parameters of i/ and ly, determined from 1991 to 1994. The 
narrow indicates the time at which the usual glitch of the 
magnitude Av jv = 8 x 10~^° was observed. This moment 
divides the observations on two parts. After that time we 
observe three slow glitches. 

During the interval from 1985 to 1994 we observed 
the distinct linear decreasing of the first derivative of 
the period. The value of the second derivative is i> = 
9 X 10-25s-3. The data from 1995 to 2006 show three 
jumps of values of v and Aj/, which can be called as 
slow. The characteristic property of slow glitches is the 
gradual increasing of the rotation frequency during 200- 
300 days without subsequent relaxation. The increasing 
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Fig. 1. The time dependance (a) of the first derivative 
of the rotation frequency v and (b) frequency resiuals Ai^ 
for the pulsar B1822— 09. Three slow glitches are observed 
from 1995 to 2005. 



of the relative rotation frequency /S.v/v due to the slow 
glitch is of 10~*. The total increasing for three glitches is 
A^/j/w 7- 10-**. 

We see that three slow glitches observed from PSR 
B1822-09 shows the quite different properties than well 
known glitches from another pulsars, for example, from 
Crab and Vela pulsars. They arc connected with processes 
occurred in a core and a crust of a neutron star. But the 
reason of the slow glitch also can't be connected with a 
magnctosphcric effects. If a pulsar magnetosphere is re- 
sponsible for a slow glitch then it must change its moment 
of inertia J on the value of A J/ J = —^vjv w — 10~®. The 
pulsar magnetosphere does not possess the such moment 
of inertia. The estimation of the maximum value of J,„ 
gives J,„ < / (i?^/87rc^)r^(ir. Here B is the magnitude of 
the magnetic field in a magnetosphere, c is the speed of 
light, r is the distance from the rotation axis. The integra- 
tion is over a pulsar's magnetosphere from a star's surface 
to the light surface = ^Ticjv. As a result we obtain 
J™ < BlR^l2^ « 5 X W^iBo/W'^GfiR/lQ^cmfgcm^. 
R is the star's radius. Thus, Jm/ J < 5 • 10"-^'^, where 
J « W^^gcm^ is the standard value of the moment inertia 
of a neutron star. So, the changes of the rotation frequency 
on the value of Aiy/iy w 10~^ as observed can be only by 
changes in a star's body. 



3. Connection between the inclination angle and 
slow glitches 

3.1. Reconstruction of a star's shape 

The neutron star shape has no the spherical symmetry. 
Due to the fast rotation its extension across the rotation 
axis is larger than that along. As a result the tensor of iner- 
tia of a star Ja^ has no equaled diagonal components. The 
value of Jzz (the axis z is directed along the rotation axis) 
is little larger than the components which are perpendic- 
ular to the rotation axis J±±: Jzz = J, J±± = J ~ SJ, 
5J<^J . The ratio 5J/ J is of the order of the ratio of the 
centrifugal force to the gravitation force 



5J_ 
1 



Gpn 



Here f2 is the frequency of a star's rotation, G is the 
gravitation constant, pn is a star's density. Substitut- 
ing the standard value of a neutron star density p = 
5 X 10^^ gcm~^ , we obtain 



6J_ 
J 



^ w 1.2 X lO^^p-^ f — 



If the inclination angle changes then the rotation axis 
deflects from the symmetry axis z (the magnetic field 
is frozen to the star's body). The axis of rotation be- 
gins to precess around the symmetry axis. That is de- 
scribed by relations fJ^^ = f2 sin a cos $; fij, = sin a sin 
Viz = ricosa = const{t); d^/dt ~ ilcosa{5J/ J), 
from which the period of precession Ppr follows Ppr = 
P{J/SJ)/cosa w 0.83 • 10^ {pn/ p) {cos a)-^s. For the 
pulsar B1822-09, which has the period P = 0.769s, the 
precession period is equal to (a<Cl): Ppr ~ 3.8 x lO'^.s = 
4.4 days. This time is significantly less than the pulsar's 
age P/P = 1.5 X lO^'^s = 5 x lO^years and the obser- 
vation time. Averaging over the precession period Ppr, 
we get < ft >= VLzez. The observed period is equal to 
P = 2tt/VIz. 

The torque acting on star results the spin down of a 
star (P > 0) and also changing of the inclination angle x- 
As a result the angle a increases. The star's shape, initially 
matched with the rotation axis, becomes asymmetric with 
respect to the instant axis of rotation. The star, rotating 
around the axis not coinciding with the axis of symmetry, 
undergoes the stress. This stress intends to reconstruct the 
star's shape in order to match it with the rotation axis. At 
a 7^ the part of the centrifugal force is not compensating 
by the gravitation force. So. there exists the force directing 
tangent to the surface of constant density (see figure 2). 

The force density /, acting on the unit volume, is 
/ = pn^^ra, where r is the distance from the axis of 
rotation. The arising tangent stress equals r = F/s ~ 

h 

an'^J pn{h'){R - h')dh' (here R is the star's radius). The 



stress T will intend to dcfiect the star matter column on 
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Fig. 2. The star rotates around the axis not coinciding 
with the symmetry axis. The angle between axises is a. 
There arises the tangent stress r intending to reconstruct 
the star's shape. 



the angle 



in- 



Pn{h') 

E{h') 



{R - h')dh'. 



The quantity E is the elastic module, it depends on 
the matter density p„, E = Es{pn/ Ps)'^^^ ■ Quantities ps 
and Eg are values of the density and the elastic module 
on the star's surface. The neutron star crust near the sur- 
face consists of the crystallic iron. Its density on the star 
surface is of ps ~ lO^g/cm^. The density increases with 
the depth /i, p„ = + [h/HY]. The quantity H is the 
characteristic scale of a crust density growth, H k, 10m 
(Blandford, Applegate, Hernquist, 1983). The value of the 
elastic module Es can be estimated from the magnitude 
of E for a steel on the Earth, i?o « 2 x lO^^din/cm^. Then 
Es = Eo{ps/pof/^ « 3 • lO^dm/cTO^ (at po « IQg/cm^). 
Resulting we obtain the following expression for the angle 
of deflection of star matter from the vertical 



m 



Es 



[R - h')dh' 

[l + (/l7H)3]l/3- 



At h>H the angle of deformation approaches the max- 
imum value oi 4) ^ a{^^ PsRH / Eg), which is much less 



crust (j{pn). The elastic limit a is proportional to the value 
of E and is of the small part of it, a{pn) = (Js^Pn/ PsY^^ 
(fT « 10"^£' on the Earth at the room temperature). The 
value CTs is the elastic limit for the neutron star crust at 
its surface. The condition t = cr defines the limiting value 
of the angle ai at which the crust matter will move under 
the action of stress. After that the star shape becomes 
symmetric with respect the rotation axis. 



(l + x3)4/3 



n^PsRH 



J{1 - Px'){l + x'^)dx' 



where x ^ h/H, (3 = H/R « IQ-^^i. The right hand 
side of this expression has the minimum at x = 2^^/"^. This 
minimum does define the limiting value of the angle ctcr- 
At this value of a the star will reconstruct its shape, and 
the matter motion will be mainly at the depth h w Q.8H, 



1.9 



n^PsRH 



(1) 



The value of the plastic limit a depends strongly on 
the matter temperature. More exactly, on the ratio of the 
binding energy of atoms to the temperature. This ratio 
is proportional to the quantity pn^^^/T. The matter den- 
sity on the star's surface ps « lO^g/cm^ is much greater 
than that for metals on the Earth po w lOg/cm^. The 
star surface temperature is also much higher, T « 10^ K. 

1 /3 

So, the parameter ps jT is of the order less than that 

1/3 

for the Earth conditions, p^ ITq{Tq « 300 if), namely 
pl'^^jT K. 6.5 X 10-2(pyVro)(T/105i^)-i. Such value of 
the binding parameter corresponds to the strongly heating 
body, for which the plastic limit a is at least of the order 
less than that for the room temperature. Because of that 
we adopt the following value for (Ts, ~ 10^^ din / crri^ . 
The obtained estimation of the critical value of the angle 
between the axis of rotation and the symmetry axis for the 
pulsar B1822-09 is acr « 2 x 10^^. The angle a can not 
exceed this limit. During the evolution of the inclination 
angle x the angle a increases from zero to Then, the 
crust reconstructs, and the angle a returns to zero value. 
This reconstruction is observed as a slow glitch. 

3. 2. Change of the star rotation frequency after the glitch 

Let us discuss the behavior the star rotation frequency 
after the glitch. The value of H. after the reconstruction 
Jl' can be determined form the condition of angular mo- 
mentum conservation Jap^fS = const. After the glitch the 
angular velocity is directed along the symmetry axis, i.e. 



than the angle a, « 3 x lO^^a. Thus, we see that the Japflp' = Jfla'- Thus, i^x' = — (5J/ J) sina cos 

fly' = — 5 J/ J) sin a sin <i>; fi^' = ilcosa. Here the 
quantity CI is the star rotation frequency before the glitch. 
Let remind that the mean rotation frequency before the 
glitch is < >= 271 /P = 57 cos a. Calculating of 17', we 



elastic deformation is not enough to reconstruct the shape 
of a star to the equilibrium form. 

While the angle a grows the tangent stress r increases 
and can become greater than the plastic limit for the iron 
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find the change of the star rotation Ail = il'— < f2 >= 
ll(aV2)(l - 25J/J) « ll(aV2). So, 

n 2 

Indeed, the observed jumps Afl for PSR B1822-09 are 
positive and are of the order of Afi/ri « 10^^. These 
values are consistent with the estimation acr ~ 10^"* , ob- 
tained above. 

3. 3. The rate of the inclination angle change 

To achieve the value a ~ 10^** — 21". and also the change 
of X on the same value it is necessary large enough time 
t « (P/ P)a. For a usual pulsar this time is of 3 x Ky^years. 
PSR B1822— 09 possesses the much higher value of the 
spin down P = 52.36 x 10~^^, and the corresponding time 
diminishes to i ~ 50years. That is also large for the ob- 
servations. But the pulsar B1822— 09 is on their stage of 
evolution that this time becomes shorter and is of sev- 
eral years. The matter is while approaching the value of 
inclination angle x to the limiting value (0 or 7r/2 de- 
pending on the mechanism of losses), when the energy 
losses become minimal, the change of x becomes larger 
than that of fl. For the magnetodipole losses the invari- 
ant is the quantity ilcosx = const. It gives the rate of Xj 
X = —[P/P) coix- At the angle x closed to zero, cotx3>l 
and x^ The time of changing of x on the value of Ax 
equals t = {P/ P)Axtanx^{P/ P)'^X- The same relation 
takes place also for the current losses. The invariant value 
for this case is ilsinx = const (x """Z^)- The angle x 
approaches 7r/2. The time is t ~ (P/P)Axcotx, and it is 
also smaller than (P/P) Ax at x ~ tt/2. 

Existence of the interpulse in radio emission of PSR 
B1822-09 tell us that the value of x is close to or to 
7r/2. Rankin (1990) determined this angle, x ~ 86°. The 
conclusion that the angle x is close to the value of 7r/2 
is confirmed by the fact that the shift between the main 
pulse and the interpulse is equal to 180° and does not 
depends on the radio frequency (Gil et al., 1994). 

One more argument in favor of orthogonal rotator 
(x ~ 7i'/2) is the measuring braking index n = flfl/fl^ 
has the large value of n « 145. Such high value of the 
braking index, n^l, means that the rate of change of the 
inclination angle x is much larger than that of the rota- 
tion frequency O. Taking into account only the rotation 
frequency evolution it is impossible to get such high value 
of 71, n « 3. In the relation (l = f{n,x) at n^l the depen- 
dence on X becomes determinant f] « xdf /dx- By this the 
braking index is n « x(^/^)(c'ln//9x)- Thus, the brak- 
ing index is defined by x ~ n{ri/Q){d\ii f /dx)~^^^/^ 
at n^l. 

For the current losses which give the evolution of 
the inclination angle to 7r/2, /(x) oc cos^x + <^^/4 (Bc- 
skin, Gurevich, Istomin, 1993). Here the quantity S is 
the angular size of the polar cap, S = (1.95Pri/c)^/'^ = 



2.3 • 10-2(i?/10fcm)i/2 w 1.3°. The value of x has 
the sharp maximum near cosx ~ S/2, at which x = 
— {(l/Q) cosx/ (cos^ X + '5^/4). The characteristic value of 
X in the interval of angles tt/2—x ~ (5 is x = —S~'^{Cl/il) = 
S~^P/P « AiP /P. At such high rate the change of x on 
the value of Ax ~ 10^"* is during the time scale of several 
years. We suggest that the pulsar B1822-09 demonstrates 
such kind of evolution. 

3.4- Comparison with observations 

We consider that observed slow glitches of B1822-09 are 
reconstructions of the neutron star shape. Wc observed 
after 1995 three slow glitches of the pulsar rotation, 
Avi/v = 1.2 X 10-*, /^v2/v = 2.85 x lO^®, Ai/g/t/ = 
2.75 X 10^^. Knowing the relation ~ 2IS.V jv^ we obtain 
the following estimations for the angles ol\ ~ 1.5 x 10^^ (« 
30"), ai = 2.4 X 10-4 (« 50"), = 2.3 x lO"* (« 48"). 
We see that observed values of a are close to the critical 
one acr = 2 x 10"^, derived above from the condition of 
plastic limit. To achieve these angles we need the time of 
t = a/x. 

For the orthogonal rotator (x ~ t^I^ the quantity x 
is X = 6~^P/P « 2.9 X lO-^^fi-^, and the time to achieve 
a = 2 X 10"^ is i = 7 X lO^s = %QMays. This time is in 
well agreement with observed intervals t2 ~ t^ k YOOdays. 

Thus, we come to the conclusion that slow glitches 
observed from PSR B1822— 09 is explained by the model 
of the neutron star reconstruction, when the asymmetry 
of the star's shape is restored by the plastic deformation 
in the crust. We also showed the pulsar B1822-09 is the 
orthogonal rotator, and x ~ 90°. This is the confirmation 
of the current losses conception. 



4. Conclusion 

We showed that the probable reason of slow glitches of the 
pulsar B1822— 09 is the fast change of its inclination angle 
X- During slow glitches the inclination angle changes on 
the value of 2' for 10 years. Significant increasing of the 
rate x takes place when the angle x approaches 90°, and 
the condition tt/2 ~ x^S is fulfilled {6 is the angular size 
of the polar cap 5 ~ {2Rfl/cy^^). In this region of angles 
X = S'^P/P-»P/P. 

For the fast evolving pulsars like Crab, Vela or 
B1509— 58, for which the ratio P-15/P is large quantity 
of 10'^ — lO''^, the time of accumulation of angle acr is 
of the order of one year. They have to demonstrate slow 
glitches. But they have came to the state when x = 90°. 
For another group of pulsars with the ratio P_i5/P«10^ 
the accumulation time is of 30 years. And only for the 
pulsar B1822-09, which has the angle x close to 90°, but 
not equal that value, the accumulation time is of several 
years. 
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